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By mapping the signal-recycling (SR) optical configuration to a three-mirror cavity, and then 
to a single detuned cavity, we express SR optomechanical dynamics, input-output relation and 
noise spectral density in terms of only three characteristic parameters: the (free) optical resonant 
frequency and decay time of the entire interferometer, and the laser power circulating in arm cav- 
ities. These parameters, and therefore the properties of the interferometer, are invariant under 
an appropriate scaling of SR-mirror reflectivity, SR detuning, arm-cavity storage time and input 
power at beamsplitter. Moreover, so far the quantum-mechanical description of laser-interferometer 
gravitational-wave detectors, including radiation-pressure effects, was only obtained at linear order 
in the transmissivity of arm-cavity internal mirrors. We relax this assumption and discuss how the 
noise spectral densities change. 
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I. INTRODUCTION 



A network of broadband ground-based laser interferometers, aimed at detecting gravitational waves (GWs) in 
the frequency band 10 — 10 4 Hz, is already operating. This network is composed of GEO, the Laser Interferometer 
Gravitational- wave Observatory (LIGO), TAMA and VIRGO (whose operation will begin in 2004) (fl. The LIGO 
Scientific Collaboration (LSC) Q is currently planning an upgrade of LIGO starting from 2008, called advanced LIGO 
or LIGO-II. Besides the improvement of the seismic isolation and suspension systems, and the increase (decrease) of 
light power (shot noise) circulating in arm cavities, the LIGO community has planned to introduce an extra mirror, 
called a signal- recycling mirror (SRM) j|, [|, at the dark-port output (see Fig. |l|). The optical system composed 
of SR cavity and arm cavities forms a composite resonant cavity, whose eigenfrequencies and quality factors can be 
controlled by the position and reflectivity of the SR mirror. These eigenfrequencies (resonances) can be exploited to 
reshape the noise curves, enabling the interferometer to work either in broadband or in narrowband configurations, 
and improving in this way the observation of specific GW astrophysical sources || . 

The initial theoretical analyses || [| and experiments |(| of SR interferometers refer to configurations with low 
laser power, for which the radiation pressure on the arm-cavity mirrors is negligible and the quantum-noise spectra 
are dominated by shot noise. When the laser power is increased, the shot noise decreases while the effect of radiation- 
pressure fluctuation increases. LIGO-II has been planned to work at a laser power for which the two effects are 
comparable in the observational band 40-200 Hz [0. Thus, to correctly describe the quantum optical noise in LIGO- 
II, the results have been complemented by a thorough investigation of the influence of radiation-pressure force on 
mirror motion [Q, ||, || |l0| . The analyses revealed that SR interferometers behave as an "optical spring" . The dynamics 
of the whole optomechanical system, composed of arm-cavity mirrors and optical field, resembles that of a free test 
mass (mirror motion) connected to a massive spring (optical fields). When the test mass and the spring are not 
connected (e.g., for very low laser power) they have their own eigenmodes: the uniform translation mode for the free 
mode and the longitudinal-wave mode for the spring. However, for LIGO-II laser power the test mass is connected 
to the massive spring and the two free modes get shifted in frequency, so the entire coupled system can resonate 
at two pairs of finite frequencies. Near these resonances the noise curve can beat the free mass standard quantum 
limit (SQL) for GW detectors |nj. Indeed, the SQL is not by itself an absolute limit, it depends on the dynamical 
properties of the test object (or probe) which we monitor. This phenomenon is not unique to SR interferometers; it 
is a generic feature of detuned cavities [|l2| [l3| and was used by Braginsky, Khalili and colleagues in conceiving the 
"optical bar" GW detectors juj. However, because the optomechanical system is by itself dynamically unstable, and 
a careful and precise study of the control system should be carried out |lCj|. 

The quantum mechanical analysis of SR interferometers given in Refs. |8j, |9[ |l0| , was built on results obtained by 
Kimble, Levin, Matsko, Thorne and Vyatchanin (KLMTV) for conventional interferometers, i.e. without SRM. 
For this reason, both the SR input-output relation ^] and the SR optomechanical dynamics JTc| were expressed in 
terms of parameters characterizing conventional interferometers, such as the storage time in the arm cavities, instead 
of parameters characterizing SR interferometers as a whole, such as the resonant frequencies and the storage time 
of the entire interferometer. Therefore, the analysis given in Refs. Ji|, |], [l0| is not fully suitable for highlighting the 
physics in SR interferometers. 

In this paper, we first map the SR interferometer into a three-mirror cavity, as originally done by Mizuno [ fl5| , though 
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FIG. 1: We draw a signal- (and power-) recycled LIGO interferometer. The laser light enters the interferometer from the left 
(bright port), through the power-recycling mirror (PRM), and get split by a 50/50 beamsplitter (BS) into the two identical (in 
absence of gravitational waves) arm cavities. Each of the arm cavities is formed by the internal test-mass mirror (ITM) and the 
end test-mass mirror (ETM). No light leaves the interferometer from below the BS (dark port), except the lights induced by 
the antisymmetric motion of the test-mass mirrors, e.g., due to a passing-by gravitational wave, or due to vacuum fluctuations 
that originally enter the interferometer from the dark port. A SRM is placed at the dark port, forming a SR cavity (marked 
by thick dashed lines) with the ITMs. 



in the low power limit and neglecting radiation-pressure effects, and by Rachmanov ]iq | in classical regimes. Then, as 
first suggested by Mizuno |15 , we regard the very short SR cavity (formed by SRM and ITM) as one (effective) mirror 
and we express input-output relation and noise spectral density ||, and optomechanical dynamics || as well, in terms 
of three characteristic parameters that have more direct physical meaning: the free optical resonant frequency and 
decay time of the entire SR interferometer, and the laser power circulating in arm cavities. By free optical resonant 
frequency and decay time we mean the real and inverse imaginary part of the (complex) optical resonant frequency 
when all the test-mass mirrors are held fixed. These parameters can then be represented in terms of the more practical 
parameters: the power transmissivity of ITM, the amplitude reflectivity of SRM, SR detuning and the input power. 
An appropriate scaling of the practical parameters can leave the characteristic parameters invariant. 

In addition, in investigating SR interferometers || ||, [To) the authors restricted the analyses to linear order in the 
transmissivity of arm-cavity internal mirrors, as also done by KLMTV |?J for conventional interferometers. In this 
paper we relax this assumption and discuss how results change. 

The outline of this paper is as follows. In Sec. |l| we explicitly work out the mapping between a SR interferometer and 
a three-mirror cavity, expressing the free oscillation frequency, decay time and laser power circulating in arm cavity, 
i.e. the characteristic parameters, in terms of SR-mirror reflectivity, SR detuning and arm-cavity storage time, which 
are the parameters used in the origin al description j^, |9| . An interesting scaling law among the practical parameters is 
then obtained. In Sees. HI and IV A the input-output relations, noi se spe ctral density and optomechanical dynamics 
are expressed in terms of those characteristic parameters, 
detuned cavity of the kind analyzed by Khalili [0| 
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In Sec. IV B we map the SR interferometer to a single 
Sec. IV C we show that correlations between shot noise 
and radiation-pressure noise in SR interferometers are equivalent to a chang e of the optomechanical dynamics, as 
discussed in a more general context by Syrtsev and Khalili In Sec. [VP , using fluctuation-dissipation theorem, 
we explain why optical spring detectors have very low intrinsic noise, and are then preferable to mechanical springs 
in measuring very tiny forces. In Sec. ^ we derive the input-output relation of SR interferometers at all orders in the 
transmissivity of internal test-mass mirrors. Finally, Sec. [v] summarizes our main conclusions. Appendix [A| contains 
definitions and notations, Appendix ^ discusses the Stokes relations in our optical system and in Appendix |c| we 
give the input-output relation including also next-to-leading order terms in the transmissivity of arm-cavity internal 
mirrors. 

In this manuscript we shall be concerned only with quantum noise, though in realistic interferometers seismic and 
thermal noises are also present. Moreover, we shall neglect optical losses [see Ref. || where optical losses in SR 
interferometers were discussed]. 
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II. DERIVATION OF SCALING LAW 
A. Equivalent three-mirror— cavity description of signal-recycled interferometer 

In Fig. [j], we draw a signal- and power- recycled LIGO interferometer. The Michelson-type optical configuration 
makes it natural to decompose the optical fields and the mechanical motion of the mirrors into modes that are either 
symmetric (i.e. equal amplitude) or antisymmetric (i.e. equal in magnitude but opposite in signs) in the two arms, 
as done in Refs. B |9|, and briefly explained in the following. In order to understand this decomposition more 
easily, let us for the moment ignore the power-recycling mirror (PRM) and the signal-recycling mirror (SRM). 

First, let us suppose all mirrors are held fixed in their equilibrium positions. The laser light, which enters the 
interferometer from the left of the beamsplitter (BS), excites stationary, monochromatic carrier light inside the two 
identical arm cavities with equal amplitudes (marked with two + signs in Fig. |l|) and thereby drives the symmetric 
mode. To maximize the carrier amplitude inside the arm cavities, the arm lengths are chosen to be on resonance 
with the laser frequency. When the carrier lights leave the two arms and recombine at the BS, they have the same 
magnitude and sign, and, as a consequence, leak out the interferometer only from the left port of the BS. No carrier 
light leaks out from the port below the BS. For this reason, the left port is called the bright port, and the port below 
the BS is called the dark port. Obviously, were there any other light that enters the bright port, it would only drive 
the symmetric mode, which would then leak out only from the bright port. Similarly, lights that enter from the dark 
port would only drive the antisymmetric optical mode, which have opposite signs at the BS (marked in Fig. Q) and 
would leak out the interferometer only from the dark port. 

Now suppose the mirrors (ITMs and ETMs) move in an antisymmetric (mechanical) mode (shown by arrows in 
Fig. [j]) such that the two arm lengths change in opposite directions — for example driven by a gravitational wave. 
This kind of motion would pump the (symmetric) carriers in the two arms into sideband lights with opposite signs, 
which lie in the antisymmetric mode, and would leak out the interferometer from the dark port (and thus can be 
detected). On the contrary, symmetric mirror motions that change the two arm lengths in the same way would induce 
sidebands in the symmetric mode, which would leave the interferometer from the bright port. Moreover, sideband 
lights inside the arm cavities, combined with the strong carrier lights, exert forces on the test masses. Since the 
carrier lights in the two arms are symmetric, sidebands in the symmetric (antisymmetric) optical mode drive only 
the symmetric (antisymmetric) mechanical modes. In this way, we have two effectively decoupled systems in our 
interferometer: (i) ingoing and outgoing bright-port optical fields, symmetric optical and mechanical modes, and (ii) 
ingoing and outgoing dark-port optical fields, antisymmetric optical and mechanical modes. 

When the PRM and SRM are present, since each of them only affects one of the bright /dark ports, the decoupling 
between the symmetric and antisymmetric modes is still valid. Nevertheless, the behavior of each of the subsystems 
becomes richer. The PRM, along with the two ITMs, forms a power recycling cavity (for symmetric optical modes, 
shown by solid lines in Fig. [l]). In practice, in order to increase the carrier amplitude inside the arm cavities H, this 
cavity is always set to be on resonance with the input laser light. More specifically, if the input laser power at the 
PRM is /in, then the power input at the BS is Iq = 4/i n /T p , and the circulating power inside the arms is I c — 2Io/T, 
where T p and T are the power transmissivities of the PRM and the ITM. The SRM, along with the two ITMs, forms 
a SR cavity (for the antisymmetric optical modes, shown by dashed lines in Fig. |l|). By adjusting the length and 
finesse of this cavity, we can modify the resonant frequency and storage time of the antisymmetric optical mode [|J , 
and affect the optomechanical dynamics of the entire interferometer [[l0| . These changes will reshape the noise curves 
of SR interferometers, and can allow them to beat the SQL ||, Qj. 

Henceforth, we focus on the subsystem made up of dark-port fields and antisymmetric optical and mechanical 
modes, in which the detected GW signal and quantum noises reside. In light of the above discussions, it is convenient 
to identify the two arm cavities as one effective arm cavity, and map the entire interferometer to a three-mirror cavity, 
as shown in Fig. || In particular, the SR cavity, formed by the SRM and ITMs is mapped into a two-mirror cavity 
(inside the dashed box of Fig. ^) or one effective ITM. The antisymmetric mechanical motions of the two real arm 
cavities is equal or opposite in sign to those of this system. The input and output fields at the dark port corresponds 
to those of the three- mirror cavity, a and b (shown in Fig. |J). Because of the presence of the BS in real interferometer 
(and the absence in effective one), the optical fields inside the two real arms is ±1/V2 times the fields in the effective 
cavity composed of the effective ITM and ETM. As a consequence, fields in this effective cavity are \/2 times as 
sensitive to mirror motions as those in the real arms, and the effective power in the effective cavity must be 

-^arm — 2/ c . (1) 

Therefore, both the carrier amplitude and the sideband amplitude in the effective cavity are \f2 times stronger than 
the ones in each real arm. In order to have the same effects on the motion of the mirrors, we must impose the effective 
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FIG. 2: We draw the three-mirror cavity which is equivalent to a SR interferometer in describing the antisymmetric opti- 
cal/mechanical modes and dark-port optical fields. The SR cavity, which is mapped into a two-mirror cavity (in dashed box) 
can be viewed as an effective mirror, with four effective reflectivities and transmissivities, p', r' (for fields entering from the 
right side), and p, r (for fields entering from the left side). The input and output fields, a and 6, corresponds to those at the 
dark-port of the real SR interferometer. 



ETM and ITM to be twice as massive as the real ones, i.e. 

m arm = 2m . (2) 

We denote by T and R = 1 — T the power transmissivity and reflectivity of the ITMs, L — 4 km is the arm length, 
and we assume the ETMs to be perfectly reflecting. The arm length is on resonance with the carrier frequency 
ujq = 1.8 x 10 15 sec -1 , i.e. ujqL/c — Nir, with TV an integer. We denote by p and I the reflectivity of the SRM and the 
length of the SR cavity, and <ft = [^oV c ]mod27r the phase gained by lights with carrier frequency upon one trip across 
the SR cavity. We assume the SR cavity to be very short (~ 10 m) compared with the arm-cavity length. Thus, we 
disregard the phase gained by lights with sideband frequency while traveling across the SR cavity, i.e. fll/c — » 0. The 
three-mirror cavity system can be broken into two parts. The effective arm cavity, which is the region to the right of 
the SR cavity, including the ETM (but excluding the ITM), where the light interacts with the mechanical motion of 
the ETM. This region is completely characterized by the circulating power I c , the arm length L and the mirror mass 
to. The (very short) SR cavity, made up of the SRM and the ITM, which does not move. This part is characterized 
by T, p and <j). 

Henceforth, we assume the radiation pressure forces acting on the ETM and ITM to be equal, and the contribution 
of the radiation-pressure-induced motion of the two mirrors to the output light, or the radiation-pressure noises due to 
the two mirrors, to be equal. [These assumptions introduce errors on the order of max{f2_L/c, T}.] As a consequence, 
we can equivalently hold the ITM fixed and assume the ETM has a reduced mass of 

,„„ = V„. (3) 



B. The scaling law in generic form 

As first noticed by Mizuno p5| , when the SR cavity is very short, we can describe it as a single effective mirror with 
frequency-independent (but complex) effective transmissivities and reflectivities (see Fig. ^) p , f (for fields entering 
from the left) and p' , r' (for fields entering from the right), and write the following equations for the annihilation (and 
creation, by taking Hermitian conjugates) operators of the electric field [see Appendix [A] for notations and definitions]: 

j±(0.)=tfk±(n)+Ta±(Q,), 6±(0) =r'fe±(0)+pa ± (0). (4) 

Among these four complex coefficients, p' , the effective reflectivity from inside the arms, determines the (free) optical 
resonant frequency loq + Q of the system through the relation: 

[Note that the carrier frequency ujq is assumed to be on resonance in the arm cavity, i.e. tooL/irc = integer.] It turns 
out that if we keep fixed the arm-cavity circulating power I c , the mirror mass to and the arm-cavity length L, the 
input-output relation (a — b) of the two-port system (|j) is completely determined by p' alone or equivalently by the 
(complex) free optical resonant frequency CI. To show this, we first redefine the ingoing and outgoing dark-port fields 




FIG. 3: We plot p and cj> - tt/2 versus e for A = 2tt x 100 Hz (solid line), 2tt x 200 Hz (dotted line), 2tv x 500 Hz (dashed line) 
and 2tt x 1000 Hz (dashed-dotted line), having fixed T = 0.033. 



as: 

a±(n) = ^-a±(n), 6±(n) = -^6±(n). (6) 
m m 

This redefinition is always possible since we can freely choose another (common) reference point for the input and 
output fields. Secondly, using the Stokes relations given in the Appendix [b|, we derive the following equations: 

J± (ft) = p' k ± (Q) + \t 1 2± (fi) = p' fc± (ft) + V1-|P'I 2 «± (") » (7) 
&±(ft) = \t\ k ± (fi) - T S± (ft) = \/l - |p' | 2 fc± (ft) - p'* a± (ft) , (8) 

from which we infer that the output fields 6±(ft) depend only on p' or equivalently on fi. Thus, if we vary the 
interferometer characteristic parameters T, p and 4> such that p 7 is preserved, the input-output relation do not 
change. We refer to the transformation among the interferometer parameters having this property as the scaling law. 



C. The scaling law in terms of interferometer parameters 

In this section we give the explicit expression of the scaling law in terms of the practical parameters of the SR 
interferometer. We start by deriving the effective transmissivities and reflectivities p, r, p' and f 1 in terms of T, 
R = 1 — T, p and <fi. By imposing transmission and reflection conditions at the ITM and SRM, and propagating the 
fields between these mirrors (see Fig. |J), we get the following equations: 

Ta±(n)+pe l,p y ± (n) = x±(Vl) , VT k±{n) - VRe l,p x±(Sl) = y±(Cl) , (9) 
- pa±(Q) + t y±(tl) =6±(ft), VRk±(fl) + y/T e*+ x±(Vl) = j±(Sl) , (10) 



where the reflection and transmission coefficients of ITM and SRM are chosen to be real, with signs {+vT, — vR}, 
{+t, — p} for light that impinges on a mirror from outside the SR cavity; and {+VT, {+t, +p} for light that 

impinges on a mirror from inside the SR cavity. Solving Eq. (||) for x± and y± in terms of a± and b± , plugging these 
expressions into Eq. ( ^L0| ) and comparing with Eq. (Q) we obtain: 



P = 1= i P = 1= i T = T = — . (11) 

l + VRpe 2 ^ l + VRpe 2 ^ l + VRpe 2 ^ 



It can be easily verified that these coefficients satisfy the Stokes relations (JB9|)— (B1C ) . The scaling law can be obtained 
by imposing that p 1 does not vary. This gives: 

' -const.. (12) 
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FIG. 4: We plot p and <J> — ty/2 versus T for three sets of optical resonances: (A, e) = (2% x 194.48 Hz, 2n x 25.42 Hz) (solid 
lines), (A,e) = (2tt x 228.10Hz, 2tt x 69.13Hz) (dotted lines) and (A,e) = (2tt x 900Hz,2tt x 30 Hz) (dashed-dotted lines). We 
mark with a square and a triangle the special configurations selected in Refs. |, |, |hj, with (T,p,4>) = (0.033, 0.9, tt/2- 0.47), 
and the current LIGO-II reference design Q, with (T, p, <f>) = (0.005, 0.96, tt/2 - 0.06), respectively. 



Using Eq. (^), we derive the (complex) free optical resonant frequency in terms of T, p and </>: 

~ ic VR + pe 2 ^ . 

il = — log == = — A — u, 13 

where we trade O for two real numbers, the resonant frequency A and decay rate (inverse decay time) e. For any 
choice of T, the parameters p and (j> can be expressed in terms of A and e by solving Eq. ([n]) in terms of pe 2l f The 
result is: 



e -2eL/c e 2i\L/c _ ^j} 
1 - ^fR e -2£L/c e 2i\L/c 



fle 2 '^- — V1L (14) 



In Fig. |H we plot p (left panel) and 4> — tt/2 (right panel) as functions of e for four typical values of A: 2ir x 100 Hz 
(solid lines), 2tt x 200 Hz (dotted lines), 2ir x 500 Hz (dashed lines) and 2ir x 1000 Hz (dashed-dotted lines), while 
fixing T — 0.033. In Fig. |[ we plot p and 4> — tt/2 as functions of T, as obtained from Eq. (|l4[), for three sets 
of optical resonances: (A, e) = (2ir x 194.5Hz, 2ir x 25.4Hz), plotted in solid lines, which goes through the point 
(T, p, (ft) = (0.033, 0.9, tt/2 — 0.47) (marked by a square), which is the configuration selected in Refs. |], 0, [n|; (A, e) = 
(2tt x 228.1 Hz, 2ir x 69.1 Hz), plotted in dotted lines, which goes through the point (T, p, (f>) = (0.005, 0.96, tt/2 - 0.06) 
(marked by a triangle), which is the current LIGO-II reference design and (A, e) = (2ir x 900 Hz, 2ir x 30 Hz), 
plotted in dashed-dotted lines, which is an example of a configuration with narrowband sensitivity around a high 
frequency. As T, p and cf> vary along these curves, the input-output relation is preserved. 

As done in Refs. 0,0], we now expand all the quantities in T and keep only the first nontrivial order. [The accuracy 
of this procedure will be discussed in Sec. For the crucial quantity p' a straightforward calculations gives: 

ff = l-- P * . (15) 
So the scaling law at linear leading order in T is: 

T TT7^ = const -- (16) 



Moreover, applying Eq. (15) to Eq. (g), we derive the following expression for the (free) optical resonant frequency at 
leading order in T: 

1 l-pe 2 ^ Tc _ -2psin20-z(l-p 2 ) 

il + pe^AL 1 + p 2 + 2pcos20 7 ' 1 ' 
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where 7 = Tc/AL is the half-bandwidth of the arm cavity. The frequency coincides with the frequency f2_ introduced 
in Ref. JhJ. [Since the authors of Ref. jl0| used the quadrature formalism, they had t o intr oduce another (free) optical 
resonant frequency which they denoted by f2+ = —171. See discussion around Eq. (A12) in Appendix [A].] Thus, at 
linear order in T we have: 



A = 



2p7 sin 2(f) 
1 + p 2 + 2p cos 2<p 



(1-P 2 )7 



1 + p 2 + 2p cos 2<p 



Finally, using Eqs. (B9) and Eq. (15) we obtain the coefficients redefining the fields a±(Q) and 6±(f2) in Eq. 

r (1 + p) cos <j) + i(l — p) sin <j> 



sj\ + 2p cos 24 + p 2 



(18) 



(19) 



III. INPUT-OUTPUT RELATION AND NOISE SPECTRAL DENSITY IN TERMS OF 

CHARACTERISTIC PARAMETERS 

A. Input— output relation 

In this section we shall express the input-output relation of SR interferometer (at leading order in T) only in terms 
of the (free) optical resonant frequency, Q = —A — i e, and the parameter l c , defined by 



8w Ic 
m Lc 

where the circulating power I c is related to the input power at BS Iq by: 



Ic - j^Io 



Using Eq. ( |l9| ) and the results derived in Appendix [A] [see Eqs. (A8), (AlO) and (All)] we transform Eqs. 
are given in terms of annihilation and creation operators, into equations for quadrature fields: 



and 



"1 

«2 



(l + p)cos(f> — (1 — p)sin</> 
V'l + 2pcos2< ? (> + p 2 V (1-P)sin0 (l + p)cos</> 



(l + p)cos0 (1 — p)sin<^ 



V / l + 2pcos20 + / 9 2 V -(1-P)sin0 (l + p)cos<£ 



ax 
a-2 



(20) 

(21) 
which 

(22) 
(23) 



Inserting the above expressions into Eqs. (2.20)-(2.24) of Ref. ||, and using Eqs. (|l8|)-(pT|), we get the input-output 
relation depending only on the characteristic or scaling invariant quantities A, e and t c : 



°11 °12 
^21 L/ 22 
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5(1) 
^2 



l SQL 



where we define: 



and 



MW = [A 2 - (n + iej 2 } SI 2 — A 
tfff = =r! 2 (fi 2 -A 2 + e 2 ) + A lc , 



C 



(i) 
12 



2eAft 2 . 



= 2e A Q 2 - 2e i c , 



(24) 

(25) 
(26) 



and 



5 W = -2A 1 /e7^0, =2(e-if2)nVeZ;, 



^SQL = 



M 2 ! 2 



(27) 



(28) 
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is the free-mass SQL for the gravitational strain h(Q) in LIGO detectors JO]. The quantity l c has the dimension of a 
frequency to the third power (O 3 ). Since it is proportional to the laser power circulating in the arm cavity, it provides 
a measure of radiation-pressure strength. In order that radiation pressure influence interferometer dynamics in the 
frequency range interesting for LIGO, we need 

L C >n 3 GW => j c > mLc "K (29) 

OUJQ 

which gives I c > 100 kW for typical LIGO-II-parameters and Oqw = 2-7T x 100 Hz. The input -output relation ( |24| ) 
is more explicit in representing interferometer properties than that given in Ref. B, and can be quite useful in the 



process of optimizing the SR optical configuration 19 . From the last term of Eq. ( j24| ) we observe that as long as 
the SR oscillation frequency A ^ 0, both quadrature fields contain the GW signal. Moreover, the resonant structure, 
discussed in Ref. |§, |[, is readily displayed in the denominator of Eq. (|24|), given by Eq. (P5|). As we shall see in 
Sec. |y|, the shot noise and radiation-pressure noise, and the fact they are correlated, can also be easily worked out 
from Eq. (H). 

In Ref. we found that one of the SR resonant frequencies, obtained by imposing M" = 0, has always a positive 
imaginary part, corresponding to an instability. This instability has an origin similar to the dynamical instability 
induced in a detuned Fabry-Perot cavity by the radiation-pressure force acting on the mirrors jl^, To suppress it, 
we proposed |l(J a feed-back control system that does not compromise the GW interferometer sensitivity. Although 
the model we used to describe the servo system may be realistic for an all-optical control loop, this might not be the 
case if an electronic servo system is implemented. However, results obtained in Refs. |2(| would suggest it does. In 
any case, a more thorough studying should be pursued to fully clarify this issue. In this paper, we always assume 
that an appropriate control system of the kind proposed in Ref. [Tcfl is used. 

Finally, when A = (which corresponds to either p = 0, or p ^ 0, <fr — 0, 7r/2) Eq. ( |24| ) simplifies to 

which exactly coincides with Eq. (16) of Ref. [Q for a conventional interferometer, but where 

(5), tf-g^ftpj- (3D 

The simple relations (3^), ([5l]) nicely unify the SR optical configuration 4> = 0, 7r/2 (denoted by ESR/ERSE in Ref. |J) 
with the conventiona -interferometer optical configuration. 

B. Noise spectral density 

The noise spectral density can be calculated as follows 0, ^] . Assuming that the quadrature — b\ sin C + &2 cos £ 
is measured, and using Eq. (0), we can express the interferometer noise as an equivalent GW Fourier component: 



h n = h SQL Ab c , (32) 



where 



~ (Cff sin C + Cff cos C) ai + (eff sin C + C$ cos Q a 2 

Ab C = ■ v33j 

D[ 1} sm( + D { 2 1} cosC 

Then the (single-sided) spectral density Sj^f), with / = S1/27T, associated with the noise h n can be computed by the 
formula [see Eq. (22) of Ref. @]: 

2Tr6(n~n')S C h (f) = (m\h n (n) hl(n') + hl(n') h n (Q)\m) . (34) 
Assuming that the input of the whole SR interferometer is in its vacuum state, i.e. |in) = |0a), and using 

(0 5 |ai(fi)at(O0 +51(00^(0)105) = 2ir 8(Sl - Sf) S tj , (35) 
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we find that Eq. (|34|) can be recast in the simple form (note that Cffi € 3?): 



(dQ sinC + C { 2 \ ] cosC) 2 + (C$ sinC + C$ cosC 



S" h = hi QL + — 3 (36) 

£>i J sin C + D 2 > cos C 

Plugging into the above expression Eqs. (|26j), fl27| ) we get the very explicit (and very simple!) expression for the noise 
spectral density: 



Q 2 h 2 9 

Si "•ROT, r,„ » \ n on r , „ . , 1 m * 



- A(e 2 + A 2 + 2e 2 cos2C)-e(e 2 -A 2 )sin2C] +^ [2e 2 (l + cos2C)-2eAsin2C + A 2 ] I. (37) 



IV. OPTOMECHANICAL DYNAMICS IN TERMS OF CHARACTERISTIC PARAMETERS 

The scaling laws ([hJ), (16) could have been equivalently derived by imposing the invariance of the optomechanical 
dynamics El . In this section we express all the relevant quantities characterizing the SR optomechanical dynamics 
in terms of the scaling invariant parameters A, e and t c . 

A. Radiation-pressure force 

In Ref. || we assumed that SR interferometers can be artificially divided into two linearly coupled, but otherwise 
independent subsystems: the probe V, which is subject to the external classical GW force G and the detector T>, 
which yields a classical output Z. The Hamiltonian of the overall system is given by [see Sec. IIB in Ref. |J for 
notations and definitions]: 

H = H v + H v -x(F + G), (38) 

where x is the operator describing the antisymmetric mode of motion of four arm-cavity mirrors and F is the radiation- 
pressure or back-action force the detector applies on the probe. In the Heisenberg picture, using the superscript (1) 
for operators evolving under the total Hamiltonian i?, and superscript (0) for operators evolving under the free 
Hamiltonian of the detector H-p, the equations of motion in Fourier domain read 0: 

Z (1 \n) = Z^(n)+R ZF (Q)x^(Q,), (39) 

f«(si) = F (0 \n) + R FF (n)x (1 Hn), (40) 
x w (n) = l h(n) + R xx {n) f w (n) . (41) 

where R XX (Q.) = — 4/m/Sl 2 , h(fl) is the gravitational strain [see Eq. (2.15) of Ref. |flo|| l related to the GW force in 
Fourier domain by G(Q) — — (to / '4) L SI 2 h( fl) , while the various Fourier-domain susceptibilities are defined by: 

f + OO 



■ r+oo 

RAB(n) = j;J dre^[A(0),B(-T)}, 



(42) 



where [A(i), B{t')\ is the commutator between operators A and B. As discussed in Sec. I, LIGO-II has been planned 
to work at a laser power for which shot noise and radiation-pressure noise are comparable in the observational band 
40-200 Hz. In Sec. IIIA of Ref. @ the radiation -pressure force was explicitly derived. Here, we want to express it, 
and the other crucial quantities entering the equations of motion (|39|)--(|4l|) in terms of the characteristic parameters 
A, e and 

l c = mi c = — . (43) 

Ij c 

Using Eqs. (|l8"|) a straightforward calculation gives the rather simple expressions: 



(0) j el c ti (in-e)a 1 (fl) + \a 2 (n) , 

[ j ~ V 2 (SI - A + ie) (n+ A + ie) ' [ ' 
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(A 2 



^ 2 )ai(fi) + 2Xea 2 (n) 



(fi - A + ie) (fl + A + ie) 
-2Aeoi(fl) + (A 2 - e 2 - fi 2 )5 2 (fi) 
(0 - A + ie) (0 + A + ie) 



(45) 
(46) 



X 



2% (SI- \ + ie){n + \ + ie) ' 

' el c (e — ifi) 

~2ft (fi - A + ie) (O + A + ie) 



(47) 
(48) 



The optical pumping field in detuned Fabry-Perot resonator converts the free test mass into an optical spring having 



very low intrinsic noise |14 
in SR interferometers, is a. 
expression is given by: 



]. The ponderomotive rigidity ifpond.rig., which characterizes the optomechanical dynamics 
so responsible of the beating of the free mass SQL [see Sec. IIIC of Ref. [fill] and its explicit 



K 



pond 



(fi) = -r ff (Q) = —Z 



A 



4 (fi - A + ie) (fi + A + ie) 



(49) 



As long as the free optical resonant frequency A differs from zero, -Kpond is always non- vanishing. Moreover, in order 
to have a (nearly) real mechanical resonant frequency at low frequency, we require A < [as can be obtained by 
imposing K pond (fl = 0) > 0.] 



B. Equivalence between noise correlations and change of dynamics 

As derived in Ref. || [h]], the output of SR interferometers, when the first or second quadrature of the outgoing 
dark-port field is measured, can also be written as: 

0i(fi) =Zi(Q) + R xx (fl) [^(fi)+G(fi)] , i = l,2 (50) 

where: 

z -w = €m- m) = F '''' in) - R ^€§)' ! ='' 2 ' (51) 

Expressing these quantities in scaling- invariant form [here the first or second quadrature refers to b\ or b 2 , so the Z\^ 2 
discussed here are related to those in Ref. by the rotation (p3|)], we get: 

■Zi(n) = J^r \ [{X 2 -e 2 -Q 2 )a 1 {Q)+2eXa 2 (^)] , (52) 

V ex c A 

Zm = Xl^r^x [2eAai(ft)-(A 2 -e 2 -Q 2 )a 2 (fi)] , (53) 

V el c (e — ill) 1 J 

and 



^(Oi = J^ M n)> (54) 



F 2 (Q) = y |^ j- [2ea!(n) - AH 2 (0))] . (55) 

The form of Eq. (|5^) , along with the fact that the operators Zi (f2) and Ti (fi) are proportional to 1/ ^[T c and \[T C , made 
it natural to refer to them [ |To| as effective output fluctuation and effective radiation-pressure force. The quantum 
noise embodied in Z^ifl) is the shot noise, while the quantum noise described by iFiifl) is the radiation-pressure or 
back-action noise. The operators Zi(Q), Ti(£l) satisfy the following commutation relations J9|, |l(j, |ll|: 

[z i (n),z}(n')] = o = [F i (ii),F}(Q.')}, [Zi(n),F}(n')] = -2Trih8(n-n'), 1 = 1,2. (56) 
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If the output quadrature i is measured, the noise spectral density (|36[), written in terms of the operators Zi and T%, 
reads 

S h ,i(fl) = jp {S Zt zM + 2n xx {n)K[Sr t zM\+K x {ti) Sr t rM} , (57) 
where the (one-sided) cross spectral density of two operators is expressible, by analogy with Eq. (|34|), as 

27r<$(n-fi') s AB (n) = (Oa\A(n)B\n') + B^{n')A(n)\o d ). (58) 

In Eq. (p7|) , the terms containing SziZi, SfiTi an d ^[^JFi-zJ should be identified as shot noise, radiation-pressure 
noise and a term proportional to the correlation between the two noises, respectively JTl[ . The noise spectral densities 
expressed in terms of the scaling invariant quantities A, e and X c are rather simple and read: 

s ZlZ An) = ^ n + X)2 + ^- X)2 + ^ , (59) 

s f01 _ 2^ [(ft + A) 2 + e 2 ] [(ft-A) 2 + £ 2 ] 

(0) - 5: 7w+m ' ( } 



S FlFl (tl) = ^, (61) 
_ K c (4 e 2 + A 2 ) 

s> a *,(«) - g7 £2 + fi2 , (62) 



SWO) = fi ^ ^ " 2) , (63) 
o ™ * A(A 2 +3e 2 -fl 2 ) 

Note that in our case Syr. Zi is real, thus Sjr. z . = Sz^- ^ i s straightforward to check that the following relation is 
also satisfied: 

Sz iZi (ft) (ft) - S Zt r t (ft) (ft) - Ti 2 , i = l,2. (65) 

Since in SR interferometers Sz^ 7^ 0, the noise spectral density Sh i is not limited by the free-mass SQL for GW 
interferometers (S'sql = ^1ql)j as derived and discussed in Refs. || El 

We want to show now that cross correlations between shot noise and radiation-pressure noise are equivalent to some 
modification of the optomechanical dynamics of the system composed of probe and detector, as originally pointed out 
by Syrtsev and Khalili in Sec. Ill of Ref. More specifically, we shall show that for linear quantum measurement 
devices, at the cost of modifying the optomechanical dynamics, the measurement process can be described in terms 
of new operators Z' and T' with zero cross correlation. 

In Ref. || the authors found that the most generic transformation which preserves the commutation relations ( p36| ) 
is of the form [see Eq. (2.25) in Ref. §]: 



Z'M\_ ia f L U L 12 \ ( Ziifi) 

- e 1 L 2i l 22 j l ^(ft) 



(66) 



with a, Lij S and detLy = 1. Under this transformation the output ( |50| ) becomes: 

Oi(Q) = e' la [L 22 - R XX (Q) L 21 ] Z&Q) + e- la [-L 12 + R xx (fl) L n ] ^(fi) + R xx (fl) G(ft) . (67) 

By imposing that the system responds in the same way to electromagnetic and gravitational forces, .F'(ft) and G(ft), 
we find the two conditions: e za = ±1 and R XX (D,) (Ln =p 1) = L\ 2 . The transformation we have to apply so that the 
correlations between new fields Z[(£i) and J 7 / (ft) are zero, give the following set of equations: 



T ( S Zi zM Sz^M \ T t _ ( s 2lzl (Q) 



(68) 
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When Szr — S?z S Si, as it happens in SR interferometers, the above conditions can be solved in infinite ways. 
A simple solution, suggested by Syrtsev and Khalili |l7|j , is obtained by taking a — and L n = 1. In this case, a 
straightforward calculation gives: L i2 = 0, L 2 i — —Sz^/ 'SziZi and L 2 2 = 1- The output becomes: 



o' i (n)=z<(n)+xfm [rm+G(m , o:-(n)=o i (o)%H, (69) 

Xi (0) 



where x? ff , the effective susceptibility, is given by: 



with 



Y eS (n) = Rxxjty (7Q) 

Xl [ ' l + R xx ^)Sz^MISz t zM' [ ' 

The spectral densities of the new operators Z[ and T[ are: 

Sz'.z'.m = s ZiZi m , = - fe^S ■ i=i > 2 ' ( 7l ) 

9 , - nic ^ 2 + » 2 ) frn 

"w" j - 2 [(ft + A) 2 + e 2 ][(ft-A) 2 + e 2 ] ■ 1 j 
These new operators satisfy the condition [see Eq. (^)] : 

Sz'.z'Xty S&r>(n) = h 2 , i = 1,2. (74) 



C. Equivalence to a single detuned cavity and frequency-dependent rigidity 



At the end of Sec. [I B we discussed under which assumptions radiation-pressure effects were included in the 
description of SR interferometers in Refs. || |l(J. There, the authors assumed that radiation pressure forces acting 
on ETM and ITM are equal, and disregarded ETM and ITM motions during the l ight round-trip time in arm cavities. 
In this case the ITM and SRM can be considered fixed, and as shown in Sec. [I A it is possible to map the SR optical 



configuration to a three-mirror cavity with only the ETM movable. We shall see explicitly in this section that, since 
the very short SR cavity can be regarded as a single effective mirror, we can further map the SR interferometer to 
a single-detuned cavity with only the ETM movable, which is exactly the system that Khalili discussed in Ref. [fj"3| . 
[More specifically, the single-detuned cavity has (complex) free optical resonant frequency ujq — X — ie, ETM mass 



./2 



and circulating power I a 



2I C . See Eqs. (hh, (g) and (§).] 



Marm ' 

If the output quadrature i is measured, the noise spectral density expressed in terms of the operators Z[ and T' il 
can be written as: 



d2 

XX 



V- 



>Z'.Z 



(75) 



In order to make explicit the connection with Ref. Jl3| , we evaluate the noise spectral density for xqw = L h/2 and 
we denote it by S XGW - It reads: 



1 



rarm J ' 



where as discussed above /j, a 



= m. By rewriting the generalized susceptibility into the form, 



4 - Marm Q 2 + 4*rf(f2) ' 

we introduce, as Khalili also did [O, the effective rigidity Kf s (il), defined by 



(76) 



(77) 



(78) 
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f(Hz) 



FIG. 5: Plot of yj Sh,2 (continuous lines) and A k /S , ™ 1 2 n (dashed lines) versus frequency / for T — 0.033, e = 2-zr x 25.0 Hz, 
A = 2n x 191.3 Hz, and two different values of the laser power circulating in the arm cavities: I c — 300 kW (lighter-color lines) 
and I c = 600 kW (darker-color lines). The free-mass SQL line (black straight line) is also shown for comparison. 



More explicitly, 



A — fi 2 



4 [(fi — A) 2 + e 2 ] [(fi + A) 2 + e 2 ] ' 



3e 2 + A 2 -fi 2 



4 [(n — A) 2 + e 2 ] [(n + A) 2 + e 2 ] 



(79) 
(80) 



Those expressions, in particular Eqs. (|76|), ( pd|) agree with those derived by Khalili [jljj for single detuned cavity 
[see Eqs. (19) and (21) in Ref. |p|] if we make the following identifications (this paper — > Khalili): A — > 5, e — > 7, 
2L/ arm /c = 4L/ c /c — > £ (energy stored in the single cavity), xf ff /4 — * X, and 4_K'| sff — > K. Note that in Ref. jl3) it is 
always assumed that the second quadrature is measured. 

The description of the measurement system in terms of the uncorrelated fields, Z[ and TL yields another way of 



understanding why in SR interferometers the free mass SQL, S\ 



SQL 



loses its significance. Indeed, by using 



Eq. ([74]), we get S z[z[ = h 2 /S^ 
obtain, 



Plugging this expression into Eq. (|75|), minimizing with respect to Sj^i^i, we 



h 



i + R* a {n)Kf (n)| ^(n) ' 



and the minimal noise spectral density is, 



2h 
Z 2 



xf(n) 



5: 



SQL 



(81) 



(82) 



which can be formally regarded as a non-free-mass SQL for the effective dynamics described by xt ■ To give an 
example, in Fig. || we plot the square root of the noise spectral densities and S™ 1 ^ versus frequency / having 
fixed e = 2-7T x 25.0 Hz, A = 27r x 191.3 Hz, for two different values of the laser power circulating in the arm cavities: 
I c = 300 kW and I c = 600 kW. For comparison we also plot the free-mass SQL line. As we can see from the plot, 
5^2" can g° quite below the free-mass SQL. 



The effective dynamics can be also used to optimize the performance of SR interferometers |lc 
following equation, 

corresponds to resonances produced by the effective rigidity, at which x° S ~* 00 and, using Eq. (pig), 



The roots of the 



(83) 



(84) 
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FIG. 6: Plot of the square root of the noise spectral density Sh versus frequency / for (i) triple-zero case (continuous line) 
with A = 2tt x 123.2 Hz, e = 27r x 13.8 Hz, and I c = 320 kW and (ii) three-single-zero case (dashed line) with A = 2ty x 191.3 Hz, 
e = 2ir x 25.0 Hz and I c — 590 kW. For comparison we also show the free-mass SQL line (black straight line). 



As observed by Khalili [[13|, we could expect that the more the roots of Eq. fl83| ) coincide, the more broadband the 
noise curve will be. For example, we could expect that interferometer configurations with double or triple zeros be 
optimal. However, as we shall see, those configurations are not much better than some of the three-single-zero cases. 

Assuming the second quadrature (i — 2) is observed, we obtain for the triple-zero case [see also Eqs. (29), (30) and 
(31) in Ref. (if]]: 




V281T 




^triple zero — \ Z ^ ■ (85) 



In Fig. ^ we plot the square root of the noise spectral density Sh .2 versus frequency / for the triple-zero case having 
fixed iltripiczoro = 27T x 100 Hz, i.e. the (free) oscillation frequency A = 2ir x 123.3Hz and e = 2tt x 13.8 Hz. The 
SQL line is also plotted. For comparison we also show the noise spectral density Sh,2 corresponding to a solution of 
Eq. ( pf ) with three-single zeros: A = 27r x 191.3 Hz, e = 2tt x 25.0 Hz and I c — 590 kW. As mentioned, the spectral 
density in the triple-zero case is not significantly broad-band, especially if compared with the three-single-zero case. 

This result originates from the non-universal nature of the curve Sf^ 1 . The SQL (p8|) does not change if we 
adjust (by varying the circulating power) the balance between shot noise and radiation-pressure noise and find the 
interferometer parameters whose noise curve can touch it. By contrast, the curve changes when we adjust (by 
varying the circulating power or the optical resonant frequencies) the effective shot and radiation-pressure noises, 
Sz'.z' an d Sjri t jr>. [The change of S™ 1 ^ as Ic is varied can be also seen from Fig. g.] As a consequence, the fact 
that S"™ 11 is low and broad-band for a certain configuration cannot guarantee the noise curve will also be optimal. In 
particular, in the triple-zero case, Eq. ( |83| ) already fixes all the interferometer parameters, leaving no freedom for the 
noise curve to really take advantage of the triple zeros. The fact that only a non-universal minimum noise spectral 
density exists in SR interferometers arises in part because of the double role played by the carrier light. Indeed, the 
latter provides the means for measurement, and therefore determines the balance between shot and radiation-pressure 
noises, but it also directly affects the optomechanical dynamics of the system, originating the optical-spring effect. 

Finally, Braginsky, Khalili and Volikov [^J have recently proposed a table-top quantum-measurement experiment 
to (i) investigate the ponderomotive rigidity effect present in single detuned cavity and (ii) beat the free mass SQL. 
Although the table-top experiment will concern physical parameters very different from LIGO-II, e.g., the test mass 
m ~ 2 x 10~ 2 g, L ~ 1 cm, f2 ~ 10 4 s _1 , I c ~ 1-10 W, however, because of the equivalence we have explicitly 
demonstrated between SR interferometers and single detuned cavities, the results of the table-top experiment could 
shed new light and investigate various features of SR optomechanical configurations relevant for LIGO-II. 



D. Optical spring equivalent to mechanical spring but at zero temperature 



When proposing the optical-bar GW detectors 0], Braginsky, Gorodetsky and Khalili pointed out that detuned 
optical pumping field in Fabry-Perot resonator can convert the free test mass into an optical spring having very low 
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intrinsic noise. In this section we illustrate why this happens in SR interferometers and why optical springs are indeed 
preferable to mechanical springs in measuring very tiny forces. 

The Heisenberg operator in Fourier domain x^ (Q) describing the antisymmetric mode of motion of SR interfer- 
ometer, satisfies the following equation [see Eqs. (|39|), (El]) above and also Eq. (2.20) of Ref. B]: 



(86) 



Using Eq. (EJ) we get: 



, . _ £ A 2 + (e - ifi) 2 

The noise spectral density associated with x is: 

S x (ty = \x(m 2 S F (n) , (88) 

where: 

7r*(n-n') s a (n) = (Oa\x^(n)x^\n')\o & ), n6(n-n') s F (n) = (o~ a \F^{n)F^°^(n')\o- a ) . (89) 

More explicitly, 

— [(n _ A)2 + e2] [(n + A)2 + e2] - ( 9 °) 

For the optical spring, which is made up of electromagnetic oscillators in their ground states (the vacuum state), we 
have [see e.g., Chapter 6 in Ref. @ §|]: 

s x (n)>2ft|9f[x(n)]|, (91) 

which can be regarded as a zero-temperature version of the fluctuation-dissipation theorem. For a mechanical system, 
e.g., a mechanical spring, with the same susceptibility, but in thermal equilibrium at temperature T 3> hCl/k, where 
k is the Boltzmann constant, the standard version of fluctuation-dissipation theorem says, 

S*(fi)=4^|9f[x(n)]|. (92) 

If we assume il ~ 2ir x 100 Hz, Ml/k ~ 5 x 10~ 9 Kelvin, the condition T 3> %Q,/k is always valid for any practical 
mechanical system. As a consequence, 

^mcch. spring _ C-opt. spring ^3) 
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AXIX 




A/271 (Hz) 2000 A72jr(Hz) 2000 

FIG. 8: We plot the fractional error AA/A (in left panel) and Ae/e (in right panel) as a function of A and e. The quantities 
AA and Ae are the difference between the value of A and e obtained from the first-order-T free optical frequency (B4J) and the 
exact one (13). 



At T — 3Q0K, n/2n = 100 Hz, we get fig"**. spring _ 10 H gopt. sp ring_ ThuS! because of the very large coefficient 
kT/hfl in Eq. (|93|), fluctuating noise in an optical spring is always much smaller than in a mechanical spring! 

For SR interferometers described in this paper, the fluctuating noise S x does not saturate the inequality in Eq. (|9l]) . 
This can be inferred from Fig. |?| where we plot R = S x (f) / (2fi versus /, where S x has been obtained from 

Eqs. (|S7|), ( |88| ) and (|90|), for the following choice of the physical parameters: m = 30 kg, T = 0.033, 7 = 2n x 98.5 Hz, 
with A = 2ir x 191.3 Hz, e = 2tt x 25.0 Hz and I c = 560 kW. The minimum of R is at the frequency corresponding to 
the (free) oscillation frequency of the SR interferometer, i.e. / m ; n = X/(2tt) = 191.3 Hz. 



V. INPUT-OUTPUT RELATION AT ALL ORDERS IN TRANSMISSIVITY OF INTERNAL 

TEST-MASS MIRRORS 

To simplify the calculation and the modeling of GW interferometers, KLMTV [[?) calculated the input-output 
relation of a conventional interferometer at leading order in T and QL/c. By taking only the leading order terms in 
T, they ignored the radiation-pressure forces acting on the ITM due to the electromagnetic field present in the cavity 
made up of ITM and BS. By limiting their analysis to the leading order in flL/c, they assumed the radiation-pressure 
forces acting on the ITM and ETM are equal. In conventional interferometers, T alone determines the half-bandwidth 
7 of the arm cavities (through 7 = Tc/AL), which fully characterizes the interferometer [see Eq. (16) in Ref. j?J and 
Eqs. ©, (gl]) above]. Moreover, since f^ow is comparable to 7 and T ~ 0.005 — 0.033, the two small quantities, 
QL/c and T are on the same order, and the accuracy in expanding the input-output relation in these two parameters 
is rather under control. [Note that if 7 ~ 2ir x 100 Hz, we have T ~ 0.033.] 

In describing SR interferometers, the authors of Refs. || ^, [l0| build on the leading-order results of Ref. 0. 
However, in SR interferometers the accuracy of expanding in T can be quite obscure, because T is not the only small 
quantity characterizing SR-interferometer performances — for example the SRM transmissivity can also be a small 
quantity. Thus, to clarify the accuracy of the expansion in T, we now derive the input-output relation at all orders 
in T, and compare with the leading order result ( pif ) || ||. The calculation is much easier if we view the SR cavity 
as a single effective mirror, as done in Sec. 0. However, in doing so, we still use the assumptions mentioned at the 



beginning of this section. See also the end of Sec. II A 



A. Free optical resonant frequencies 

It is interesting to investigate the error in the prediction of the (free) optical resonant frequency introduced by using 
only the leading order terms in T and QL/c. For a generic set of T, p and </>, it can be quite complicated to characterize 
that error. For example, when p > V~R and (j> ~ 7r /2, p 1 is near —1 (in the complex plane) and the expansion jl5| ) 
around p' = 1 totally breaks down. However, we are only concerned with those parameters meaningful for a GW 
detector, and thus we limit our analysis to the region where |f2| = \/X 2 + e 2 ~ flew < 10 4 s _1 , corresponding to 
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|fii/c| ;$ 0.1. In this way \p' — 1| is always relatively small. To test the accuracy, we fix T, and for each fi 



we solve Eq. (|13|) for p and <p. Then, we insert these values into Eq. (|17|) to get the first-order-T expression for fl, 
which we denote by fjW. The result is: 



= n 




1 / flL 



o 



(94) 



From this equation we infer that since |fiL/c| < 0.1, and T is smaller than a few percents, the error in the (free) 
optical resonant frequency is not very significant (less than a few percents). In Fig. |§| we plot the fractional differences 
(denoted by AA/A and Ae/e) between the real and imaginary parts of OS 1 ' and fl, as functions of e and A for T — 0.033. 
The fractional differences are always smaller than 2.5%. 



B. Input— output relation and noise spectral density 



Using the formalism of Sec. [n| and Appendix it is rather easy to derive the exact input-output relation in terms 
of A, e and t c . The input-output relation (j-k) of the arm cavity composed of the effective ITM and ETM is: 



_ e 2iQL/c 



1 



£ ) +e mL / c y/2IC 



arm J, arm I 1 
"SQL V 



where 



_ 8/ arm cjQ _ 16/ c o;o m _ 

~ /iarrnOV _ mQV ' SQL ~ 



87/ 



/i arm f2 2 _L 2 V mfl 2 L 2 



81) 



Writing Eqs. (|7]) and (^) in terms of quadratures, that is 



and 



a i 

«2 



.., . cos ip — sin ip 
' I simp cos?/; 



cos ip sin ip 
— sin ip cos V 1 



k 2 



a i 

a 2 



(95) 



(96) 



(97) 



(98) 



where ip — arg(p'), and using Eq. (p"l|), we obtain the input-output relation (a — 6) of the three-mirror cavity, and 
thus that of the equivalent SR interferometer. They can be represented in the same form as Eq. (||), with MW, 
Cjj , and Z)| replaced by: 



M ex 
and 



4L 2 



^ _ e 2i(n+\+ie)L/c 



I _ e 2i(f2-A+ie)L/c 



Ccx 
11 — °: 

°12 

rex 
21 

5f 



22 _ 4L 2 
ft 2 c 2 



e 2i(n+A+ie)L/c _ 2i(f2- A+ie)L/c 



, (99) 



4L 2 
» 2 c 2 
4L 2 

r» 2 c 2 



1 - 2e- 2ei / c cos(2AL/c)cos(20L/c) + e~ 4£i / c cos(4AL/c) + -^ e - 4ei/c sin(4AL/c) WjOO) 

J \Vc J 

2e- 2ei / c sin(2Ai/c) [cos(2fiL/c) - e~ 2tL ' c cos(2AL/c)l + ^^e~ itL ' c sin 2 (2AL/c)l , (101) 
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2e~ inL/c ~ 2e~ 2eL/c e inL/c cos(2AL/c) 



(1 - e- 4£i / c )t CJ L 



fi 2 c 



1 - e- 4eL / c cos 2 (2AL/c)l \ , (102) 

(103) 
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FIG. 9: Comparison of first-order T-expanded (dashed line) and exact (continuous line) noise spectral density \fSu versus 
frequency /. In the left panel we use the parameters T — 0.033, p = 0.9 and <f> = 7r/2 — 0.47, m = 30 kg, and I c — 592 kW and 
show the curves for the two orthogonal quadratures bi (lighter-color lines) and 62 (darker-color lines). In the right panel we 
use T = 0.005, p = 0.964, cj> = tt/2 - 0.06, m = 40 kg, I c = 840 kW, and C = l-13ir. 



In order to compare with the results obtained in Refs. ||, ^, [T^], we have also to relate a, b to a and b. The exact 
transformations [to be compared with Eqs. (|||), (p3|)] are: 

... Am _„ .i /1 /m„-„j.\ /_ \ 

(105) 



a\ \ _ 1 / (1 + pVR) cos0 — (1 — p\/R) sin(/)\ ( a\ 



a2 J \/l + 2 P Vlcos20 + p2 i? V (l-pVl)sin^ (l + pVl)cos 



^ 7 • /l + 2 P Vi?cos20 + p2 i? V -(1 - ^) sin^ (1 + p^) cos (f> J \ b 2 



a 2 

yii-zpvficus/?)-i-p-n ... 

and 

h\ 1 ( (l + py/R)cos<f> (1 -pVR)sin<f> \f h \ ( 1Q6 ) 

As an example, we compare in the left panel of Fig. [)]the exact and first-order T-expanded noise spectral densities for 
the two orthogonal quadratures b\ and 62, having fixed T = 0.033, p = 0.9, = ir/2~ 0.47, m = 30kg and I c — 592 kW 
(which corresponds to Iq — Isql at BS) as used in Refs. j^, ^, ^0|. The T-expanded noise spectral density is given 
by Eq. (37), where we used for A, e and the redefined output quadratures Eqs. (|l8|), (|l9|). The exact noise spectral 
density is obtained from Eq. @ by replacing M«, c\]\ and t3 2 (1) with M cx , C" 3 /, and 5° x . Fr om Fig. |9[ we see 
that there is a discernible difference. In the right panel of Fig. ^, we compare the exact and first-order T-expanded 
noise spectral densities using the reference-design parameters of LIGO-II [ fl8| : T = 0.005, p = 0.964, <p = n/2 — 0.06, 
ui — 40 kg, I c = 840 kW, and Q — 1.137T. In this case, the two curves agree nicely with each other, presumably, because 
T is rather small. In the general case, if we want to trust the leading order calculation, it is not obvious how small 
T can be, since p and <fi have to change along with T to preserve the invariance of interferometer performance. For 
this reason, it is more convenient to seek an expansion that is also scaling invariant, i.e. whose accuracy only depends 

1 /3 

on the scaling-invariant properties of the interferometer. To this respect, the set of quantities \L/c, eL/c, l c L/c 
and QL/c, which are all small and on the same order, is a good choice. It is then meaningful to expand with respect 
to these quantities and take the leading order terms. We denote the noise spectral density obtained in this way by 

1 /3 

first-order X-e-L c -expanded noise spectral density. [This technique of identifying and expanding in small quantities of 
the same order can be very convenient and powerful in the analysis of complicated interferometer configurations, e.g., 
the speed meter interferometer ^2|.] Not surprisingly, doing so gives us right away the scaling-invariant input-output 

relation (p4]). In the left and right panels of Fig. |l^ we compare the exact and first-order A-e-ic^-expanded noise 
spectral densities for the two orthogonal quadratures 61^, with the same parameters used in Fig. i.e. T = 0.033, 

p = 0.9 and = tt/2 - 0.47, m = 30 kg, and I c = 592 kW (left panel) and T = 0.005, p = 0.964, <j> = tt/2 - 0.06, 

1/3 

rn = 40kg, I c = 840 kW, and £ = 1.137T (right panel). The first-order A -e-t c -expanded noise spectral density is 
obtained using for A, e and the redefined output quadratures Eqs. (|l^), (105). The agreement between the exact 



1/3 

and first-order \-e-i c -expanded noise spectral densities is much better than the agreement between the exact and 



T-expanded noise spectral densities, given in Fig. ^. 



1Q 




FIG. 10: Comparison of first-order \-e-i c -expanded (dashed line) and exact (continuous line) noise spectral density V Sh 
versus frequency /. In the left panel we use T = 0.033, p = 0.9 and <f> = n/2 — 0.47, m = 30 kg, and I c — 592 kW, and show 
the curves for the two orthogonal quadratures b\ (lighter-color lines) and i>2 (darker-color lines). In the right panel we use 
T = 0.005, p = 0.964, (j> = vr/2- 0.06, m = 40 kg, I c = 840 kW, and £ = 1.13tt. 



When either XL/c, eL/c, tc L/c or VLL/c is not small enough, the first-order X-e-L^ 3 expansion fails. An interesting 
example of astrophysical relevance is the configuration with large A and small e, which has narrowband sensitivities 
centered around a high (optical) resonant frequency. In the left panel of Fig. [n^we compare the first-order X-e-i^J 3 - 

expanded noise spectral density with the exact one, for the two quadratures 61.2 having fixed: A = 2ir x 900 Hz, 

1/3 

e = 20 Hz, m = 30 kg and I c — 600 kW. Near the lower optomechanical resonant frequency, the first-order X-e-u c 
expansion deviates from the exact one by significant amounts. However, it is sufficient to expand up to the second 
order in XL/c, eL/c, l^J 3 L/c and flL/c to get a much better agreement, as we infer from the right panel of Fig. 
[The input-output relation expanded at second order is given in Appendix 



VI. CONCLUSIONS 



In this paper we showed that, under the assumptions used to describe SR interferometers j|, ^|, [T^], i.e. radiation 
pressure forces acting on ETMs and ITMs equal, and ETM and ITM motions neglected during the light round-trip 
time in arm cavities, the SR cavity can be viewed as a single effective (fixed) mirror located at the ITM position. We 
then explicitly map the SR optical configuration to a three-mirror cavity [ll| [see e.g., Sec. |J or even a single 



detuned cavity |l3j [see Sec. IV B|. The mapping has revealed an interesting scaling law present in SR interferometers. 
By varying the SRM reflectivity p, the SR detuning <f> and the ITM transmissivity T in such a way that the circulating 
power I c and the (free) optical resonant frequency (or more specifically its real and imaginary parts A and e) remain 
fixed [see Eq. (|l8|)1, the input-output relation and the optomechanical dynamics remain invariant. 

We expressed the input-output relation (|24|), noise spectral density ( |36"| ) and all quantities characterizing the 
optomechanical dynamics, such as the radiation-pressure force ([44| ) and ponderomotive rigidity (^9|), in terms of 
the scaling invariant quantities or characteristic parameters. The various formulas are much simpler than the ones 
obtained in the original description ||, |9[ |o) . The scaling invariant formalism will be certainly useful in the process 



of optimizing the SR optical configuration of LIGO-II |19 and for investigating advanced LIGO configurations 



Moreover, the equivalance we explicitly showed between SR interferometer and single detuned cavity, could also make 
the table-top experiments of the kind recently suggested in Ref. EM more relevant to the development of LIGO-II. 

In this paper we also evaluated the input-output relation for SR interferometers at all orders in the transmissivity 
of ITMs [see Sec. 0. So far, the calculations were limited to the leading order. We found that the differences between 
leading-order and all-order noise spectral densities for broadband configurations of advanced LIGO do not differ much 
[see Fig. However for narrowband configurations, which have an astrophysical interest, the differences can be quite 
noticeable [see left panel of Fig. p!l| . In any case, we showed that by using the (very simple) next-to- leading-order 
input-output relation, explicitly derived in Appendix |^, we can recover the all-order results with very high accuracy 
[see right panel of Fig. Il[] . 

Finally, it will be rather interesting to investigate how the results change if we relax the assumption of disregarding 
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FIG. 11: For the two orthogonal quadratures fei (lighter-color lines) and 62 (darker-color lines) we compare the first-order 
A-e-tc^-expanded noise spectral density (dashed line) with the exact (continuous line) noise spectral density (in left panel) and 
the second-order A-e-z,y 3 -expanded noise spectral density (dashed line) with the exact (continuous line) noise spectral density 
(in right panel). For all the cases we fix A = 2n x 900 Hz, e = 2n x 20 Hz, m = 30kg and I c = 600 kW. 



(i) the motion of ITMs and ETMs during the light round-trip time in arm cavities and (ii) the radiation-pressure 
forces on ITMs due to light power present in the cavity composed of ITM and BS. This analysis is left for future work. 
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APPENDIX A: USEFUL RELATIONS IN THE QUADRATURE FORMALISM 

As in Refs. jj], P we describe the interferometer's light by the electric field evaluated on the optic axis, i.e. on 
the center of light beam. Correspondingly, the electric fields that we write down will be functions of time only. All 
dependence on spatial position will be suppressed from our formulae. 

The input field at the bright port of the beam splitter, which is assumed to be infinitesimally thin, is a carrier 
field, described by a coherent state with power I and (angular) frequency u)q. We denote by /gw = Q/2tt the GW 
frequency, which lies in the range 10— 10 4 Hz. The interaction of a gravitational wave with the optical system produces 
sideband frequencies luq ± Q in the electromagnetic field at the dark-port output. We describe the quantum optics 
inside the interferometer using the two-photon formalism developed by Caves and Schumaker |E3|. The quantized 
electromagnetic field in the Heisenberg picture evaluated at some fixed point on the optic axis is 



Ac Jo 2tt 

where h.c. means Hermitian conjugate and we denoted a+(f2) = a Wo+ [2 and ci-(Q) = a WQ _Q. Here A is the effective 
cross sec tional area of the laser beam and c is the speed of light. The annihilation and creation operators a±(fl) in 



Eq. (Al) satisfy the commutation relations: 

[a+,a^,] = 2tt 5(fl - Q') , [a-, at,] = 2ttS(Q - O') , (A2) 
[a + , a + >] = = [a_, a_<] , [at , a\,\ = = [al, a^,] , [a + , a_/] = = [a + , a}_,] . (A3) 
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Following the Caves-Schumaker two-photon formalism p3fl , we introduce the amplitudes of the two-photon modes as 

a + + al_ a + -a f _ 

ai = —?T' a2 = ^r ; (A4) 

a\ and a 2 are called quadrature fields and they satisfy the commutation relations: 

[01, a\,] — — [d2i Q>\i] — 2m5(Q — f2') , 

[ai,a{,] = = [ai,av] , [a 2 , 4'] = = \ a 2, a>2>] ■ (A5) 



The electric field (Al) in terms of the quadratures reads: 

E{a,i]t) = cos(wot)£i(oi;t) + sin(w t) E 2 {a 2 ;t) , (A6) 

where: 



47rftw /" + °V __ int t int . df2 



E j {a i ;t) = S j—^f J (a je -^ + a]e««)— j = l,2. (A7) 
Any linear relation among the fields a±(f2) of the kind: 

6±(fi) = /±(fi)a±(n), /+(n) = /(wo + n), /_(fi) = /(wo-fi), (as) 

can be transformed into the following relation among the quadrature fields: 



M - If (/+ + /-) *(/+-/*) A f <*i 

62 J - 2 I -* (/+ + /*) (/+ + /*) J [02 



(A9) 



In general, the above equation can be very complicated. In this paper we restrict ourselves to two special cases. The 
first case is when | /+ 1 = | /_ | and we write 

/±(0) = F(Q) e ^±( fi ) Vft >0, (A10) 



M - svo^(*+-*-)/ 2 ( cos ^ sin — \ I fl i 



and Eq. ( |A9[ ) becomes: 

(S J = ^ - ^ j ^ - j ■ (AH) 

It is easily checked that the input-output relation for the following processes: (i) free propagation in space, (ii) 
reflection and transmission from a thin mirror, (iii) reflection and transmission from one (or more) Fabry-Perot cavity 
for which ujq is either resonant or antircsonant, and (iv) reflection and transmission from one (or more) FP cavity 
whose bandwidth is much larger than the range of values of O we are interested in [in this case /(SI) can be considered 
as a constant (complex) number] are all special cases (or linear combinations) of the relation ( All ). 

The second case of interest for us is when there is one resonance at u>q + f2 r , with f2 r complex. In this case f(fl) is 
of the form: 

/M = gH „ , (A12) 

where g(oj) does not have poles. For f2 > 0, we have 

_ g(uj + 0) g*(q> - 0) 

(n + n;) + »)-(»- n r ) > - n) 
/+ + / - (0 - n r ) (n + n r *) ' 1 j 

(n + si*) g (u; + n) + (n - n r ) g > - n) 

Since the quadrature field at f2 mixes the frequencies Wo + ^ an d w o — the single resonant frequency r2 r appears in 
the above equation as a pair of resonant frequencies {Sl r , — O*}. 



and thus 
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FIG. 12: A two-port linear optical system can always be expressed in terms of four effective transmissivities and reflectivities, 
p' , r (for fields entering from the right side), and p, r (for fields entering from the left side). By taking the complex conjugates 
of the field amplitudes and inverting their propagation directions, a new set of fields related by the same set of transmissivities 
and reflectivities is obtained. 



APPENDIX B: THE STOKES RELATIONS 



The transmission and reflection coefficients of a system of mirrors, or more generally of a two-port linear optical 
system, can always be expressed in terms of four effective transmissivities and reflectivities: jo, r, p' and r' [see Fig. p"2"| . 
These quantities are generally frequency dependent (complex) numbers. For the fields shown in Fig. [l2], we have: 



p'ku. 



TCL U 



(Bl) 
(B2) 



Imposing that the two-port linear optical system satisfies the conservation of energy, we have: 

\p\' 



-' 2 < |fl 2 = l. 



IpT + |t'| 2 = i 



(B3) 



If we take the complex conjugates of all the complex amplitudes and revert their propagation directions, the resulting 
configuration is also a solution of the optical system, in the sense that the new fields are also related by the same 
sets of effective transmissivities and reflectivities. Thus, the system is invariant under time reversal. By applying 
explicitly this symmetry, it is straightforward to derive: 



PP* 
P P 



T T 
4- T T 



= o, 

p = 0. 



(B4) 

(B5) 
(B6) 



Equations (B3)-(B5) are the well-known Stokes relations |2J|. If we rewrite the transmissivity and reflectivity coeffi- 
cients as 



r e" 



and insert them into the Stokes relations (B4)-(B5), we obtain 



\P\ = \?\ 



\P\' 2 



e iv _ e iu' gi(/i+/i') _ _ e 2il/ 



lrl 2 = l: 



(B7) 
(B8) 

(B9) 
(BIO) 



APPENDIX C: INPUT-OUTPUT RELATIONS AT SECOND ORDER IN TRANSMISSIVITY OF 

INTERNAL TEST MASSES 

1/3 

The input-output relation expanded up to second order in XL/c, eL/c, i c L/c and VLL/c can be obtained in a 
straightforward way by expanding Eqs. (|99|)-( |l04|) . The new coefficients c\f and £>f ) are very simple. In fact, 

they can be represented in terms of the first-order ones, C\f and given by Eqs. (p5|)-(p7|), through the 

following formulas (truncated at the next-to-leading order): 



M (2) = (1 - 2eL/c)M w , 



(CI) 
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( C[? C[f \ ( 1 XL/c\ ( C$ \ ( 1 XL/c\ 

and 

= /=,(.A /C T)(1P)' 

It is quite remarkable that, at second order, the optomechanical resonances, determined by = 0, remain un- 

changed with respect to the first order result obtained imposing = 0. Apart from a (frequency-independent) 

rotation of the quadrature phases, the input-output relation at next-to-leading order are very similar to the leading- 
order one. 
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